function [Xmin Xmax Ymin Ymax IntensityImageData VelDirOneVelData BW L stats b CropRegion XTextSpacer YTextSpacer]...
    = CropAndProp(TempIntensityImageData, TempVelDirOneVelData, BW, countInt, ZoomOut, t, b)

% dont want to explain this right now


L = double(BW(:,:,1));
stats = regionprops(L,'basic');

Xmin        = stats.BoundingBox(1)-ZoomOut;
Xmax        = (stats.BoundingBox(1)+stats.BoundingBox(3))+ZoomOut;
XTextSpacer = (Xmax-Xmin)/16;

Ymin        = stats.BoundingBox(2)-ZoomOut;
Ymax        = (stats.BoundingBox(2)+stats.BoundingBox(4))+ZoomOut;
YTextSpacer = (Ymax-Ymin)/13;

% %    had to comment this out until figure out how to correctly index
% %    find the bounding box/crop region
%     Xmin = min(nonzeros(ROI(:,1,:))) - ZoomOut + 1;
%     Xmax = max(nonzeros(ROI(:,1,:))) + ZoomOut + 1;
%     Ymin = min(nonzeros(ROI(:,2,:,ROInumber))) - ZoomOut + 1;
%     Ymax = max(nonzeros(ROI(:,2,:,ROInumber))) + ZoomOut + 2;
%     YTextSpacer = (Ymax-Ymin)/13;
%     XTextSpacer = (Xmax-Xmin)/16;

% crop the image and reset the bounding values, xmin, max...
CropRegion         = [Xmin Ymin (Xmax-Xmin)-1 (Ymax-Ymin)-1];
TempIntensityImageData = imcrop(TempIntensityImageData, CropRegion);
TempVelDirOneVelData   = imcrop(TempVelDirOneVelData, CropRegion);

[Rsize Csize] = size(TempIntensityImageData);% Preallocate Matrices (2 speed during load?)
IntensityImageData(:,:,:) = zeros(Rsize,Csize,countInt);
VelDirOneVelData(:,:,:)   = zeros(Rsize,Csize,countInt);
IntensityImageData(:,:,1) = TempIntensityImageData;
VelDirOneVelData(:,:,1)   = TempVelDirOneVelData;

L = double(BW);  % eventually want to make this bwlabel to work with n ROIs

for i = 1:t

    BWtemp(:,:,i)     = imcrop(BW(:,:,i), CropRegion);
    Ltemp(:,:,i)      = imcrop(L(:,:,i), CropRegion);

    b(:,1,i) = b(:,1,i) - Xmin + 0.5 ; % 0.5 spacer is to place line in center of pixel
    b(:,2,i) = b(:,2,i) - Ymin + 0.5 ;

end

clear BW L
BW = BWtemp;
L  = Ltemp;
clear BWtemp Ltemp;

for i = 1:t
    stats(i) = regionprops(L(:,:,i),'basic');
end

Xmax        = Xmax - Xmin;
Xmin        = 1;
Ymax        = Ymax - Ymin;
Ymin        = 1;

end
